Dynamic percolation and Slow Relaxation in Glass-like Materials 
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Glass-like materials are nonequilibrium systems where the relaxation time may exceed reason- 
able time scales of observations. In the present paper a dynamic percolation model is introduced in 
order to explain the principal properties of glass-like materials near the dynamic transition. Here, 
contrary to conventional percolation problems, clusters are groups of particles dynamically corre- 
lated. Introducing a size dependent relaxation time and the scaling hypothesis for the distribution 
of dynamically correlated clusters the two step relaxation predicted by the mode coupling theory 
and observed in experiments is obtained. 
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I. INTRODUCTION 



In common second order transitions the new phase has 
as distinguish property the existence of long range order, 
leading to macroscopic structures like ferromagnetic do- 
mains. However, while it is believed that glass-like tran- 
sitions are second order transitions, the glass-like state 
does not exhibit the conventional long range order but 
a short range order, at a mesoscopic scale, is preferred. 
The relaxation time of these systems becomes as long 
as, or longer than, the observation time and, therefore, 
the system cannot reach equilibrium within observation 
time. Thus, it is not clear if the glass-like transition 
is a true thermodynamic phase transition or a slow re- 
laxation phenomcHa, usually referred to as dynamic or 
freezing transitionEI. 

A major step toward resolving the dilemma would be 
the determination of an order parameter underlying the 
transition. The discovery of a diverging length associated 
with the glass-like transition, similar perhaps to the cor- 
relation length in a para-ferromagnetic transition would 
enhance our understanding of glass-like transitions and 
of the glassy state. From the structural point of view the 
glassy state is characterized by some local (quenched) dis- 
order and frustrated bonds which inhibit the long range 
order and lead to the formation of clusters of mesoscopic 
sizes. Thus, any theory which try to explain the principal 
features of glass-like materials must take into account the 
existence of such clusters- In this direction several mod- 
els have been proposedH~Q. 

During the past few years a new theoretical ap- 
proach, the mode coupling theory (MCT), has been 
developed wijiph has motivated new experimental 
investigational 113. This theory provides a generalized 
kinetic equation approach to the density fluctuation dy- 
namics in supercooled fluids. It predicts two step re- 
laxation processes (a and /3 relaxations). The a relax- 
ation process exhibits a time scale which diverges at a 
crossover temperature Tc, which is somewhat above the 
calorimetric glass transition temperature Tg, and struc- 
tural relaxation experiences rapid slowing down. On the 



other hand the time scale of the j3 process is considerably 
shorter, does not diverge at Tc, and the P process can 
therefore also be observed below Tc- These predictions 
are directly accessible to experimental verificatiiC 
have been tested for many glassforming materialstlrll^ 

Recently we have presented a fractal model which ex- 
plains |the hysteresis loss effect observed in supercooled 
liquidsQ. The determinant influence of the fractal prop- 
erties of the amorphous structure in the system dynam- 
ics was emphasized introducing a fractal distribution of 
cluster sizes. However, in that occasion, the nature and 
origin of this fractal clusters was not analyzed. 

In the present paper a cluster model for the relaxation 
in glass-like materials is proposed following the general 
ideas presented inQ. As in previous works the relaxation 
of the system is obtained by superposition of indepen- 
dent relaxing clusters with a size dependent relaxation 
time. Nevertheless, in our model, differing from other 
formulations, a cluster exists through a series of relax- 
ation process giving rise to a dynamically correlated clus- 
ter (DCC). In this way the slow relaxation dynamics and 
other features of glass-like materials are explained. Re- 
sults are compared with experimental data reported in 
the literature. 

The paper is organized as follows. In section 2 the for- 
mation of DCC is analyzed. The main ideas and concepts 
about DCC and cluster dynamics are introduced. Then, 
the influence of the distribution of DCC on relaxation 
dynamics is investigated. It is obtained that the relax- 
ation follows a two step process which may be identified 
with the a and /? processes. Comparison with experi- 
mental data and MCT predictions is presented in section 
3. Finally the conclusions are given in section 4. 



II. DYNAMIC PERCOLATION MODEL 



A. DCC relaxation time 



Molecular dynamics simulations of soft-sphere 
mixture£3 have revealed that in addition to the stochas- 
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tic jump motion of single atoms correlated jump motions, 
in which a cluster of several atoms jump at successive 
closed times by permuting their positions, take place. 
Similar results have been obtained in molecular dynamics 
simulations of water, a frustrated system with multiple 
random hydrogen bond network struat4|ii:eEt3. Besides, 
MonteCarlo simulations of spin glasse^^fia also revealed 
the existence of such cooperative relaxation dynamics. 

These results suggest that the cooperative motion of 
particles (atoms, spins, etc.) in glass-like materials does 
not take place as the motion of frozen clusters but as suc- 
cessive concatenations of single particle motions. Hence, 
in this sense, a cluster is formed through the perturba- 
tion propagation from one particle to its neighbors. This 
kind of cluster will be referred here as a dynamically cor- 
related cluster (DCC). For short times only small clusters 
participate in the relaxation dynamics however, with in- 
creasing time, larger and larger clusters are incorporated. 

Let us denote a characteristic physical quantity of the 
«-th DCC particle by Xi. It may be the displacements 
around an equilibrium position in ordinary glasses, the 
magnetic moment in a spin glass, or other. Xi will be 
random variable following certain distribution, with zero 
mean value and variance (xf) = S^, 6 being a typical 
value of fluctuations. Moreover, let tq be a characteristic 
time scale of the fluctuations in Xi. 

The different DCC configurations will be then charac- 
terized by the sum 



(1) 



where s is the number of particles in the cluster, 1 > 
s < oo. In general the fluctuations Xi are correlated 
but, in a first approximation, they may be considered in- 
dependcfU;. In such an approximation the central limit 
theoremE3 tells us that, for large s, Ps{x) the limit dis- 
tribution of the sum in equation follows a Gaussian 
distribution with variance 5s^/^, i.e. 
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(2) 



Ps (x) represents the relative number of microstates or 
configurations associated to a value of x, for a cluster of 
size s. The entropy of this 'temporally isolated' subsys- 
tem is thus given by 

k 

Ss(x) = klnPs(x) = ---T^r^^ + const. , (3) 
2(3 -^s 



and the force associated to an entropy change 



h{x)=T 



dSsjx) 
dx 
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Now, from the kinetic point of view, the relaxation of 
the DCC may be reduced to the Brownian motion of a 
harmonic oscillator with generalized coordinate x, un- 
der the elastic force in equation (0) and with diffusion 



coefficient D = jr^. In this niean field approach the 
normalized relaxation is given byEII 
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with the size dependent relaxation time 
kTl 



-JT- = TQS 

U c 



(5) 



(6) 



Thus, in this mean field approximation the relaxation 
time of a DCC is proportional to its size and a pertur- 
bation of the equilibrium state relaxes exponentially to- 
wards equilibrium. As will be seen, the existence of a 
broad distribution modifies the relaxation rate. 



B. Relaxation dynamics 

Let Us be the cluster numbers, i.e. the number of DCC 
with size s per unit lattice. We assume that, in spite of 
the dynamic nature of the clusters, the scaling assump- 
tion about the cluster number still holds. Therefore, if 
the system is close to percolation e^{l — T/Tc) ~ the 
scaling assumption establishes thatca 



(7) 



where t and a are scaling exponents and f{x) is a cutoff 
fmiction (/(z < 1) « 1, and f{z > 1) < 1). Thus, the 
averaged normalized decay will be given by 



(8) 



Close to percolation the sum is dominated by the con- 
tribution of large cluster and we may substitute the sum 
by an integral. Thus, substituing equations (H) and 
(0) in equation with the variable change z = t/Tg, it 
is obtained 
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(9) 



where = To|e|^-'^/'^ and a — t ~ 2 

The first thing we note from this expression is the exis- 
tence of two characteristic times tq and . The former is 
a characteristic time for the microscopic motion while the 
last one is the characteristic time of the largest clusters 



with size 



-l/a 



It is tentative to associate these characteristic times 
with the (3 and a relaxation processes, respectively. 
Moreover, the distribution of cluster sizes is divided in 
two regions. The initial part which follows a power law 
decay and a second part with a sharp cutoff. The cutoff 
size is of the order of e"^/'^ with characteristic time r^. 
Hence, divides the time scale in two regions where the 
relaxation function will show different behaviors. These 
two time scales are analyzed below. 



2 



1. p relaxation 



with 



For t <^ Te equation (g) is approximated by 



(10) 



where 7(0;, a) is the incomplete Gamma Euler function. 
Moreover, for t <C tq, using the series expansion of 7(x, a) 
for smaU x, we obtain 
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while for t ^ To, ^{I/tq, a) ~ r(a), resulting 



0(t)«r(l + a)(- 
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(11) 
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Thus, the initial decay, due to single particle motions, is 
exponential. However, for t ^ tq cooperative dynamics 
leads to the formation of DCC, with power law distribu- 
tion sizes, giving the power tail in equation (|l2|). The 
fractal behavior in the distribution of DCC sizes leads to 
fractal properties in the time decay. 

The imaginary part of the susceptibility x" is related 
to the relaxation function through the equation 



(13) 



where FT denotes the real part of the Fourier transform. 
Then, for tq <C t <C t^, using equation ( p^ ) it is obtained 



a sec 



(14) 



This frequency dependence corresponds with the (3 relax- 
ation. The characteristic time for these time scale is tq 
which does not depend on e. Hence, if we assume that 
the divergences near Tc are given by the percolative na- 
ture of the DCC, To will not diverge at T = Tc as it is 
expected for the /? process. 



2. a relaxation 

Now let us investigate the influence of the cutoff in the 
distribution of cluster sizes in the long time (t ^ t^) re- 
laxation dynamics. In this case the form of the cutoff 
function f{x) is determinant and cannot be ruUed out. 
We will assume that this cutoff function is of the form 



/(a;) = exp(-a;^±) 



(15) 



where 0_ and 9+ stand for e < and e > 0, respectively. 

An asymptotic expansion, for t 3> t^, of (/)(t) defined in 
equation (^) with f{x) in equation ( p^ ) give the estimate 
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P± = 



a0± 
l + a9± 



(17) 



(18) 



Ca.e± is a constant which does not depend on e and, there- 
fore, Te± diverges as Te near Tc- 

Equation (|lq ) is the well known Kolraush stretched ex- 
ponential often used to fit the experimental data in the 
a relaxation process. Moreover, the characteristic time 
Te± oc Te diverges at T — Tc as it is expected for the a 
relaxation process. 



3. Static viscosity 

The static viscosity is related to the mean cluster size 
through the expression 

77 = G(t,> -Gto (5) , (19) 

where G is an elasticity modulus and (s) is given by 

(s) = 51 s'"-/ II sn, . (20) 
For T - Tc the mean cluster size defined- in equation 



J) will diverge according to (s) ~ |e| and, there- 



fore. 



T] ^ Te Te± 
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Hence, Te± the characteristic time for the a relaxation 
process shows the same temperature dependence of the 
static viscosity. 



III. DISCUSSION 

The present model pchibits the two step relaxation pre- 
dicted by the MCTETEI. The exponent a of the /3 process 
{x ~ w") is identified here with a = t — 2. More- 
over, since 2_< t < 5/2 as it is obtained in conventional 
percolationEa then < a < 1/2, which is in the range 
predicted by the MCT. On the other hand, the Kolraush 
stretched exponent of the a process, given by equation 
(|T8|), satisfies < /3± < 1 which is also the range pre- 
dicted by the MCT. Moreover, the characteristic time of 
the a process has the same temperature dependence as 
the static viscosity, in agreement with the MCT. 

Thus, our model seems to arrive to the same results ob- 
tained by the MCT, which have been extensively testecUn 
light scattering and neutron scattering experimentsllirEil. 
In order to perform a more accurate comparison let us 
investigate the behavior of x around the minimum be- 
tween the a peak and the (3 process. This minimum 
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should be near the frequency 1/re since it marks the 
transition from the P to the a process in our model. In 
this frequency range we will have the contribution of the 
(3 process with relax according to (t/ro)~" and of the 
a process relaxing as exp[— (t/rgi)''*] « 1 — {t/T^±Y^. 
Thus, a gross estimate of the relaxation function in the 
time scale will be given by considering both contributions 
which, resulting 



(22) 



where b — f3± — a. Substituting this expression in equa- 
tion (p^) and assuming a, 5 <C 1 it is obtained 



la/tT 



(23) 



This expression constitutes an interpolation formula of- 
ten used to fit the experimental data. The values of a 
and b obtained from the fit to experimental data are in 
general small, thus our supposition a, & ^ 1 is consistent 
with experiments. From this expression we obtain that 
the position of the minimum cUmin and the value of x {^) 
at the minimum Xmin satisfies 



\a/a 



(24) 



(25) 



lia/a — 1/2 then we obtain the MCT predic±iaii.B 
is consistent with the light scattering datall3'LjEj tiJ. 
However, under this assumption the static viscosity in 
equation ( ^ ) will diverge near Tc according to ~ 
with 7 = l/2a, while the MCT predicts 7 = l/(2a) + 
1/(26). Thus, when a/a = 1/2 we obtain the MCT pre- 
dictions, except for the exponent 7. Our result is more 
general since the ratio a/ a is not necessarily restricted to 
this value. 

The value of a obtained from the fit to light scatter- 
ing data for different glass-forming materials is found 
temperature independent and around 0.3. Therefore, 
T = 2 + a ^ 2.3 is slightly larger than the one obtained in 
percolation theory in three dimensions r ~ 2.15Eil. Be- 
sides, as was shown above, a/cr = 1/2 in order to obtain 
a good agreemen with the MCT and experiments. In 
the present model a/a = (t — '2)la is the critical expo- 
nent (3 of conventional percolatioi:ica which is also slightly 
smaller (0.4) in three dimensionsEll. These results sug- 
gest that the scaling exponent t is larger than the one 
obtained in conventional percolation theory. The scaling 
assumption holds but the power decay is stronger. 

On the other hand, the MCT predicts that the expo- 
nent P of the a-process does not depend on temp er a ty xe. 
However, it is found to be temperature dependentilj'EjO. 
It takes a constant value above Tc and then decreases 
near Tc, taking again a constant value below Tc- This 
high and low temperature limits are identified in our 



model with /3_ and respectively. The main contri- 
bution to this temperature dependence is given by 9± 
in equation-( T^) . For instance, in conventional percola- 
tion theoryElT6'_ = 1 (/?_ = 1/2) and ae+ = 1 - 1/d 
(/3+ — {d — l)/{2d ~ 1)), where d is the space dimen- 
sionality, and therefore /3- > /?+. While this behavior 
is in agreement with the experimental observations the 
values of f3±, using conventional percolation exponents, 
are smaller than those observed in experiments. Thus, we 
found again that conventional percolation theory can not 
give a precise quantitative description of the properties 
of glass-like materials. However, the scaling assumptions 
are still valid, but with different scaling exponents. To 
obtain a good agreement with the experimental observa- 
tions a0± should take larger values which implies that, in 
glass- forming materials, the cutoff for larger DCC sizes 
is stronger than the one expected in conventional perco- 
lation. 

This sharp cutoff for large sizes may be attributed to 
frustration effects which destroy the long range order. 
The existence of frustration, which is an inherent prop- 
erty of glass-like materials, may lead to different critical 
exponents. The percolation problem including frustra- 
tion effects (frustrated percolation) is a subject of recent 
study, and exhibits properties whici, are different from 
the corresponding unfrustrated caseEa. 

Finally we want to mention that other authors have 
used_the idea of DUC, for instance tha models by Domb 
et aa, Cohen et a£l, and ChamberlinO. However, these 
models does not analyze the properties of glass-like sys- 
tems around the dynamic transition, which is the main 
contribution of the present work. 



IV. CONCLUSIONS 

The complex dynamics of glass-forming materials was 
reduced to the formation of independent dynamic cor- 
related clusters (DCC) with a size dependent relaxation 
time. The relaxation of the system was obtained as a 
superposition of exponential relaxations over the distri- 
bution of the DCC sizes. In this way we have obtained 
the two step relaxation process observed in supercooled 
liquids. 

Our model contains as a particular case the MCT pre- 
dictions, it is more simple and it is based in a very general 
principle, the scaling assumption for the distribution of 
DCC near Tc- Besides, the present model gives the cor- 
rect temperature dependence for the Kolraush exponent 
(3, while the MCT predicts a constant value in disagree- 
ment with the experimental bahavior. 

The comparison of our results with the experimental 
data for glass-materials suggest that the percolation crit- 
ical exponent r should be larger than the one obtained in 
conventional percolation theory. Moreover, the cutoff for 
larger DCC sizes should also be stronger. These results 
were attributted to the existence of frustration effects. 
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which inhibit the formation of long range order dynamic 
structures. 

We conclude that the use of geometrical descriptions, 
the distribution of DCC size for instance, constitute an 
alternative approach to describe the dynamics of glass- 
like materials. The properties of these systems near the 
dynamic transition can be predicted using conventional 
tools of standard critical phenomena such as the scaling 
hypothesis. 
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